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ABSTRACT 

We investigate escape of cosmic ray (CR) electrons from a supernova remnant (SNR) 
to interstellar space. We show that CR electrons escape in order from high energies 
to low energies like CR nuclei, while the escape starts later than the beginning of the 
Sedov phase at an SNR age of 10^ — 7 x 10'^ yrs and the maximum energy of runaway CR 
electrons is below the knee about 0.3 — 50 TeV because unlike CR nuclei, CR electrons 
lose their energy due to synchrotron radiation. Highest energy CR electrons will be 
directly probed by AMS-02, CALET, CTA and LHAASO experiments, or have been 
already detected by H.E.S.S. and MAGIC as a cutoff in the CR electron spectrum. 
Furthermore, we also calculate the spatial distribution of runaway CR electrons and 
their radiation spectra around SNRs. Contrary to common belief, maximum-energy 
photons of synchrotron radiation around 1 keV are emitted by runaway CR electrons 
which have been caught up by the shock. Inverse Compton scattering by runaway 
CR electrons can dominate the gamma-ray emission from runaway CR nuclei via pion 
decay, and both are detectable by CTA and LHAASO as clues to the CR origin and the 
amplification of magnetic fiuctuations around the SNR. We also discuss middle-aged 
and/or old SNRs as unidentified very-high-energy gamma-ray sources. 

Key words: acceleration of particles - cosmic rays - gamma rays - shock waves - 
supernova remnants. 



1 INTRODUCTION 

I The origin of cosmic rays (CRs) is a longstanding problem 
in astrophysics. Supernova remnants (SNRs) are thought to 

, be the origin of Galactic CR nuclei and electrons. The most 
popular SNR acceleration mech a nism is the diffusive shock 
acceleration (Axford et al." 19771 : [Kirmskvl 1 19771 : iBelJ 1 19781 : 
[eiand ford fc Ostrikcr 1978). In fact, Fermi and AGILE show 
that middle-aged SNRs int eracting with mo l ecular clouds 
emit GeV gamma rays (e.g. lAbdo et al.ll2009l : [ Tavani et al.l 
I2OI0I ) and the origin of the GeV gamma rays can be inter- 
prete d as the decay of neutral pions produced by CR nuclei 
re.g.. lOhiraet al.l201ll ). In addition. X-ray observations pro- 
vide an evidence that electrons a re accelerated to high ly rela- 
tivistic energies in SNR shocks ama et al ][l99i). SNRs 
have been a lso observed by ground-based Cherenkov tele - 
scopes (e.g. iMuraishi et al 1 l2000l : [Ah aronian et al.ll2006al '). 
However, it is still unclear whether the origin of TeV gamma 
rays is inverse Compton scattering by CR electrons or the 
decay of neutral pions produced by inelastic collisions be- 
tween CR protons and ambient thermal nuclei. In addition, 
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there are many unidentified very-high-energy gamma-ray 
sources in our Galaxy, and their emission mechanism is also 
still unclear. 

Escape of CR nuclei fr om SNRs has been in vesti- 
gated by several authors (e.g., Ptuskin fc Zirakashvil: 
Ohira et al.ll2010l : ICaprioh et al.ll2010l : I Ohira 



loki 



2005; 



2011 



Drurvll201lh , and e mission from runaway CR nuclei has also 



been investigated (Aharonian & Atovan 'l996';'Gabici et ajj 
[2009. : .Ohira et al. 2011, 2012; E llison fc Bykov 2011). Run- 
away CR nuclei can emit gamma rays from the exterior 
of accelerators and its radiation spectrum is softer than 
that fro m the interior because of ene rgy-dependent diffusion 
of CRs (* Aharonian fc Atovan|[l99^ ). The escape process of 
CRs from accelerators is important because it changes the 
CR spectrum. The runaway CR spectrum depends not only 
on the CR spectrum at accelerators but also on the evolu- 
tion of the maximu m energy and the ev olution of the number 
of accelerated CRs jOhira et al.ll201ol '). In addition, consid- 
ering the escape effect, one can explain the spectral differ- 
ence between CR protons and helium observed by CREAM 
jAhn et al.ll201ol : lOhira fc Iokall201ll ) . 

However, the same processes for CR electrons suffering 
cooling effects have never been investigated so far. Cooling 
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processes via synchrotron radiation and inverse Compton 
scattering are important for CR electrons because their cool- 
ing times can be smaller than other characteristic timescales. 
Therefore, it is unclear whether CR electrons can escape 
from SNRs or not, and whether the electron cooling affects 
the spectrum of runaway CRs electrons or not, as in the en- 
ergy spectrum of CR electrons inside the SNRs. Runaway 
CR electrons can produce gamma rays outside accelerators. 
In the low density region, these gamma rays probably dom- 
inate over those from runaway CR nuclei and may be ob- 
served as unidentified very-high-energy gamma-ray sources. 
This is because CR nuclei can not produce sufficient gamma 
rays there electrons can always produce enough gamma rays 
via inverse Compton scattering of CMB photons (see equa- 
tion (7) of [Katz & Waxmann 2008). Furthermore, distant 
SNRs can not be easily identified by radio observations be- 
cause of brighter emission from molecular clouds than syn- 
chrotron radiation. 

Direct measurements of the CR electron spectrum 
may also bring us important information on escape of CR 
electrons. So far, Fermi, H.E.S.S. and MAGIC have re- 
vealed the spectrum of CR electr o ns and positrons up to 
a few TeV (lAharonian et al.|[2009l : lAckermann et ahllioiol : 



iTridon et al. 12011 '). Future experiments, such as AMS-02, 
CALET, CTA, and LHAASO, wiU measure the spectrum 
of CR electrons up to 1 — 100 TeV with good accuracy 
llKouninel 2010': Torii et al. 11 20081 : ICTA consortiumll201Gl : \CaA 
l2010l ).|Kobavashi ct al. (200^ pointed out that a few nearby 
sources like Vela and Cygnus loop may leave their own signa- 
tures in the TeV energy band, and that we will be able to see 
a spectral shape o f CR electrons and posi trons from a single 
source. Moreover, iKawanaka et al.l l|201ll ) have pointed out 
that escape of CR electrons can be investigated by future ob- 
servations of CR electrons, in particular, with a low-energy 
spectral cutoff. 

In this paper, we investigate escape of CR electrons 
from an SNR. Our stance in this paper is that because 
the evolution of the magnetic field has not been understood 
theoretically, we predict observable quantities by using phe- 
nomenological approaches and restrict the phenomenolog- 
ical models by comparing the predicted values and obser- 
vations. Assuming the evolution of the maximum energy of 
CR protons ( hereafter we treat nuclei as protons) d uring the 
Sedov phase l|Gabici et al.ll2009l : lOhira et al]|2010l ). we ob- 
tain the evolution of the diffusion coefficient, the magnetic 
field, and the maximum energy limited by synchrotron cool- 
ing (section [2]). Moreover, we calculate spatial distributions 
of runaway CR electrons around the SNR (section |3]), and 
radiation spectra from runaway CR electrons (section 
Finally, we discuss our results (section [5|. 



2 ESCAPE OF CR ELECTRONS 
2.1 Evolution of SNRs 

In this paper, in order to understand essential features of 
escape of CR electrons, we assume simple evolutions of the 
shock radius, Rsh{t), and the shock velocity, Ush(t), as in the 



following forms: 
Rsb{t) = Rsx \ 



Ush{t) 



Rs 

ts 



( t sS ts ) 



(1) 



(2) 



where t is the SNR age, and Rs and ts are the SNR shock ra- 
dius and the SNR age at the beginning of the Sedov phase, 
respectively (see Table [1] for summary) . The Sedov phase 
starts when the swept-up mass, 47rnisM'7ipi?g/3 becomes 
comparable with the ejecta mass Mej, where nisM and rup 
are the number density of the interstellar medium (ISM) and 
the proton mass, respectively. The free expansion velocity, 
itfrco = Rs/ts, is obtained from _Esn = AfojWfj(,c/2, where 
-EsN is the explosion energy of a supernova. Then, Rs and 
ts are represented by 



Rs 



ts 



2.13 pc 



(3) 



209 yr 



/ -EsN 



V 1051 erg 



/ ?^ISM \ 

U cm-V 
_ 1 

' I Mcj \ '' ( nisM Y^(A\ 
Vlcm-s^ ^' 



V IM, 



For simplicity, we assume i?s = 2 pc and ts ~ 200 yr in this 
paper. 

2.2 Relevant timescales 

The maximum energy of accelerated particles is limited by a 
finite SNR age, their cooling, or escape. Hence it is obtained 
by comparisons of timescales, which are given as functions 
of a CR energy, E, and the SNR age, t, (see Table [T|. The 
acceleration time of DSA, tacc(i5,t), is represented by 

D{E,t) 



t^cc{E,t) = ria. 



(5) 



Wah(t)^ 

where D{E,t) is the diffusion coefficient around the shock, 
and 77acc ~ 10 is a numerical factor which depends on 
the shoc k compressi on ratio and the spatial dependence of 
D{E,t) (|Drurvll983l ). First, we assume the Bohm-type diffu- 
sion of relativistic particles although the diffusion coefficient 
is sti ll unc lear when the magnetic field is strongly amplified 
(e.g. iRevill e ct al. 2008). The Boh m-type diffusion means 
that the mean free path of CRs is proportional to the gyro- 
radius of CRs, where the constant of proportionality, %(*), 
is the gyrofactor and rjg — 1 for the Bohm limit. Then, the 
diffusion coefficient of CRs is represented by 

cE 

where c, E, e, and B{t) are the velocity of light, the energy of 
CRs, the elementary charge and the magnetic field strength 
in the upstream region, respectively. In section 12.41 we will 
give explicit forms of time dependence of B{t) and rig{t) 
[see Eqs. (|18[) and (|19p ]. Then, the acceleration time can be 
expressed by 

cE 



tacc(i5, t) = '7acc??g(t)T 



(7) 



by 



' 3eB{t)u,i,{ty 
The escape time due to diffusion, tasc{E,t), is written 
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^csc(-£';^) — f]csc 



cE 



(8) 



DiE,t) 

where 77osc(< 1) is a numerical factor. 

The energy loss of CR protons above 1 GeV is due to 
the pion production by inelastic collisions, and the cooling 
time of CR protons, tcooi,p(i), is represented by 

1 



tc 



At) 



(9) 



where n and app « 3 x 10"^" cm^ are the number density 
(~ 4nisM in the downstream region) and the cross section 
of the nuclear interaction, respectively. A factor of 0.5 is the 
inelasticity of the nuclear interaction. Cooling of CR elec- 
trons above 1 GeV is due to synchrotron emission, inverse 
Compton scattering, or bremsstrahlung emission. For typi- 
cal SNRs, the energy loss is dominated by synchrotron emis- 
sion in the downstream region when the maximum energy is 
limited by cooling. The cooling time of CR electrons due to 
synchrotron emission in the downstream region, tcooi,c{E, t), 
is represented by 



(10) 



4e'^Bi{tyE ' 

where me is the electron mass and Bd{t) is the magnetic 
field strength in the downstream region. In this paper, we 
assume -Bd(i) ~ 4-B(t) because of the shock compression. 

2.3 Maximum energy of CR protons 

For CR protons, cooling is usually not important for a de- 
termination of the maximum energy. One can obtain the 
age- limited maximum energy, -Em,agc(i), from the condition 
tacc{E, t) — t as 



En 



3eB{t)Rl 

?7acc??g(t)cts 



( t ^ is ) 



(11) 



while the escape-limited maximum energy, Em,csc{t), is ob- 
tained from the condition tacc(i5, i) ~ tcsc(E,t) as 



En 



cTjaccEi-n^agc {t) 



(12) 



per 
is, ^/rf, 
cape 



Note that i5m,agc(t) is compara ble to -Em. es c (f) because 
-^T;Gsc??acc is on the order of unity (|Drurvll201lh . In this pa- 
we assume -Em,ago(i) = i5m,esc(i) for simplicity (that 
'esc??acc = 1)- If the maximum energy is limited by es- 
CRs above -Em, esc (t) leave the shock front upstream. 
However, they are caught up by the shock during the free 
expansion phase because the diffusion length is proportional 
to t^^^ but the shock radius is proportional to t. Therefore, 
the maximum energy is limited by the SNR age during the 
free expansion phase. On the other hand, CRs can leave 
the shock front during the Sedov phase because -Rsh oc t^^^ , 
which is slower than diffusion. Therefore, the maximum en- 
ergy is limited by escape during the Sedov phase. 

As shown in Equations (jlip and (|12p . the max- 
imum energy depends on the evolution of the mag- 
netic field around the shock (e.g. IPtsuskin fc Zirakashvilil 
I2OO3I ). Although the magnetic field amplificat ion a round 

20041: 



the shock is studied by linear analvses (iBel 



Rcvillc ct al.' '2007'; 'O hira et al] '2009b': Ohira fc T akaharal 
2010 ; Bykov et al 201 ll: ISchure fc Bell 20n) and numer- 
ical simulations ( Lucek fc Belli 2000l : iGiacalone fc Jokipiil 



l 2007| MNiemiec et al.' '2008*; 'Riquelme & Spitkovskv '2009|; 
llnoue et al. 2009j^ Ohha et al.. 2009a : . Vladim irov ct al. 201} 



Gargate et al.l i201ol ). the saturation of the magnetic-field 
amplification and the diffusion coefficient have not been un- 
derstood yet in detail. Here we use a phenomenological ap- 
proach based on the assumption that young SNRs are re- 
spon sible for observed CRs below the knee (|Gabici et al.l 
l2009h . The maximum energy of CR protons, -Em,p(t), is ex- 
pected to increase up to the knee energy 
until the beginning of the Sedov phase ts and decreases from 
that epoch. We may assume a functional form of -Em,p(t) to 
be 



Eni,p{t) — -Eknc 



(13) 



where a is a parameter to describe the evolution of the max- 
imum energy during the Sedov phase. In this paper, we as- 
sume that -Bm.p ~ 1 GeV at the end of the Sedov phase 
t — lO'^'^tscdov to reproduce Galactic CRs, then a — 2.6. 
Hereafter, we adopt a = 2.6. This assumption is consistent 
with the following fact. So far, about 300 radio SNRs have 
been observe d in a part of our Galax y | Case fc Bhattacharval 
Il998h . From[c ase fc Bhattacharval (|l998h . we expect that 
there are about 5 x 300 radio SNRs in our Galaxy. Assuming 
the supernova rate of 0.03 yr~^, the lifetime of radio SNRs 
is about 5 x 10* yr which is comparable to the end time 
of the Sedov phase. Radio synchrotron photons are emit- 
ted by electrons with energies of a few GeV, so that we can 
expect that the lifetime of radio SNRs corresponds to the 
escape time of 1 GeV particles. Therefore, the assumption. 



-Em,p(i = 10 ■ tSedo 



1 GeV, is reasonable. The tran- 



sition time from the Sedov phase to the radiat ive phase is 
somewhat ambiguous. iTruelove fc McKeej (|l999l ) and jPetruki 
(|2005h thought the end of the Sedov phase is 1.2 x 10* yr 
and 3 x 10* yr, respectively. In these cases, a becomes 3.66 
and 2.99, respectively. Future direct observations of CR elec- 
trons, AMS-02, CALET, CTA and L HAASO wiU be able to 
provide useful informations about a l|Kawanaka et al]|201ll : 
iThoudam fc Horande]||2012f ). 

Assuming Equation (|13|) . the age-limited maximum en- 
ergy -Em, ago for t !^ ts can be expressed by 



En 



— -S'kncc 



(14) 



while for t ^ ts, the escaped-limited maximum energy -Em,c 
can be expressed by 



-Em, CSC (^) — -Ekn 



(15) 



Moreover, from Equation (|15[) . the time when CRs with an 
energy of E escape from the SNR, Tcsc(E), can be expressed 
by 



E 



E'k.n 



(16) 



From Equations ([T} and (|16[) . the escape radius of CRs, that 
is, the SNR radius at which CRs with an energy of E escape, 
can be expressed by 



7?csc(S) = Rs 



E 



Ei, 



(17) 
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2.4 Magnetic field 

Equations (|12p . and (|13[) determine tlie evolution of 

B{t)/rig{t). During tiie free-expansion phase, botli tlie up- 
stream magnetic field and the gyrofactor are constant with 
time: B{t < ts) = -Bfree and rig{t < ts) = ??g,frcc, because 
these quantities may depend on the shock velocity which 
is constant in this phase. For the Sedov phase, we assume 
that the upstream magnetic field strength is B oc as 
long as it is larger than the value of the ISM, Bism, and 
that after the end time of the magnetic field amplification 
ffl at which B{tB.) = -Bism, B is equal to -Bism- Then, the 
upstream magnetic field, B(t), and the gyrofactor, rig{t), are 
given by 



B{t) 



and 



-^frcc 

Bfreo (t^) 

Bism 



(ts^t^tB) 
( tB < t ) 



B 



3eRl 



174 



IQi 



eV / V 200 yr 



2 pc 



tB = ts 



Ob 



a - 

_9_ 



( for 77g = ?7g,frGC 

( for B^ oc ) 
( for B^ cc u^^ ) 



tB = 



1.09 X 10^ yr 
1.82 X 10'' 
1.74 X 10^ 



yr 
yr 



( for 
( for 
( for 



B^ cc ult, ) 
B^ oc ) 



E m,cool 




(18) 



tc ts te tB t 



(19) 



where Bfrco is the amplified magnetic field during the free 
expansion phase and given by 

ee 



Figure 1. Schematic picture of evolutions of the maximum en- 
ergy and the break energy of CR electrons for rjg = »?g,free and 
tc < tB - The black, blue, and green lines show the age-limited, the 
cooling-limited, and the escape-limited maximum energy, respec- 
tively. The red dashed line shows the maximum energy of CR 
electrons. The cyan line shows the break energy owing to syn- 
chrotron cooling, ts, te, and te are the SNR age at the beginning 
of the Sedov phase, the start time of escape of CR electrons 

, and the end time of the magnetic field amplification, respectively. 
Em,e{te) is the maximum energy of runaway CR electrons. 



TT^ ' (20) 



and 7?g,frec ~ 1 is the gyrofactor during the free expansion 
phase, and the end time of the magnetic field amplification 
tfl is given by 



y Bism J ^ 

The downstream magnetic field is Bd — 4:B{rac = 697 
during the free expansion phase. In this paper, we consider 
the following three evolution models of B(t) for ts < t < tB, 



(22) 




The first model originates from the assumption that the gy- 
rofactor, rjg, is constant during B > Bism. In this case, 
B^ oc t~^'^ oc ufh for Q = 2.6. The second model is proposed 
20041 ) . The third model originates from the assump- 
tion that the pressure of the amplified magnet ic fie ld is pro - 
portional to the shock ram pressure (e.g. IVolk et al.ll2005h . 
Observations of the veloc i ty de pendence of the magnetic 
field can be found in IVinkI l|2008l ) that suggests B^ oc w^- 

From Equations (|20p - ((22]), the end time of the mag- 
netic field amplification, tB, are given by 



(23) 



where we assume a — 2.6 and Bism ~ 3 /iG. Note that 
we do not discuss the density dependence of the magnetic 



Figure 2. The same as Figure[T] but for oc m^j^ or w^j^. 



field in this paper. Furthermore, there is no guarantee that 
the magnetic field evolution is expressed by a pow er law 
form (|Ptsuskin fc Zirakashviiill2003l : lYan et al.ll2012l '). How- 
ever, the Galactic CR spectrum is approximated by a sin- 
gle power law below the knee energy, that is, CR observa- 
tions suggest absence of characteristic scale. Moreover, the 
power-la w behavior of the gyrofactor, -pg oc t°'~^^^ , is ex - 
pected bv lPtsuskin fc Zirakashldiil l|2003l ): lYan et all (|2012h . 
Therefore, we assume that all evolutions have power law de- 
pendences. We believe that this approximation is useful to 
extract essential features. 
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Table 1. Quantities and parameters appearing in Section 2 



Characteristic times 



tacc(S,t) 


acceleration time 


Eqs. lO, Q 


tcsc(S,t) 


escape time due to diffusion 


Eq. ^ 


^cool,p (^) 


cooling time of CR protons 


Eq. (ig 


tcool,c{-E'i t) 


cooling time of CR electrons 


Eq. JlOll 


tc 


start time of cooling 


Eq. ^ 


ts = 200 yr 


start time of Sedov phase 


Eqs. CJ, 6} 


tc 


start time of escape of CR electrons 


Eq. ^ 




end time of magnetic field amplification 


Eqs. ((2lJ, ll23ll 


Tcsc(-E) 


time when CRs with an energy of E escape from the SNR 


Eq. lig 



Characteristic energies 



-^m,age ('f ) 


given by Ucc {E,t) = t 


Eqs 


. CTJ, 


m 


£'m,csc{t) 


given by tacc{E,t) = tesc{E,t) 


Eqs 






-^m,cool (^) 


given by Ucc{E,t) = tcool(-E,t) 


Eq. 








given by tcool{E,t) = t 


Eq. 


l|33]l 




Em, pit) 


maximum energy of CR protons 


Eq. 


in 




Eui,c{t) 


maximum energy of CR electrons {mm{Em,!igc{t) , i?m,csc(t), £^m,cool(t)}) 


Eqs 




ll27)l. ll3T|l 


Em,S = -£'m,cool(t < is) 


cooling-limited maximum energy during the free expansion phase 


Eq. 


ll25ll 




Eh.S = Ey^its) 


break energy at the beginning of the Sedov phase 


Eq. 


l|3l 





Other quantities evolving with time 



-Rsh{*) 


SNR shock radius 


Eq. Q 




shock velocity 


Eq. lO 


D{E,t) 


diffusion coefficient around shocks 


Eqs. (O, lO 


Bit) 


upstream magnetic field 


Eqs. I©, |[18) 


Biit)=4B{t) 


downstream magnetic field 




Vg{t) 


gyrofactor 


Eqs. I©, ligi 



Others 



t 


SNR age 






E 


CR energy 






e 


electron charge 






c 


velocity of light 






Rs = 2 pc 


SNR radius at 4g 




Eqs. Q, id 


Resc{E) 


escape radius of CRs 




Eq. |[17} 


DismiE) 


diffusion coefficient in ISM 




Eq. glj 




maximum energy of CR protons at tg 




Eq. ^ 


BisM = 3 AiG 


ISM magnetic field strength 




Eqs. Jgrj 


Bfrec = 174 /.G 


amplified field strength during the free expansion phase 




Eqs. ((ISJ, ^ 




gyrofactor during the free expansion phase 




Eq. ((igj 


»?acc = 10 


numerical factor for the acceleration time 




Eq. lO 


7?csc =0.1 


numerical factor for the escape time 




Eq. ^ 


"ISM = 1 cm-3 


ISM density 




Eq. © 


Ob 


temporal decay index of the magnetic field (B oc t"'*^) 




Eqs. ((ISJ, 


a = 2.6 


temporal decay index of the maximum energy of CR protons (i?m,p oc t~°') 




Eqs. ([13J, IHSll 


/3 = 0.6 (t > ts) 


temporal index of the number of CRs electrons inside SNRs {dN/dE cntl^E 


-n 


Eq. OSJ 


s = 2.0 


spectral index of CR electrons inside SNRs (dN/dE oc E~") 




Eq. (Eg} 




normalization of the spectrum of runaway CR electrons (dN^sc/dE = AE~^ 




Eq. (ISSj 



2.5 Maximum energy of CR electrons 



The maximum energy of CR electrons is limited by the SNR 
age at early stages. As the maximum energy increases, syn- 
chrotron cooling becomes significant and limits the maxi- 
mum energy. From the condition ta_cc{E,t) = tcooi.c{E,t) 
and Equations p8|) and (|19|l . we obtain the cooling-limited 



maximum energy, i5m,cooi(t), as 
1 

£m,cool(t) ~ Em,S X \ ijs^ 



2C.R-Q-1 



[t ^ ts) 

(ts t ^ ts) 

''^ (tB =^ t) 



where Em,s is the cooling-limited maximum energy during 
the free expansion phase (t < ts) and given by 



En 



o ,3/2pl/2 
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2.02 X 



f— 



-Bkn 



eV 



is 



200 yr 



2^) ^2^) 



Em,s can be obtained from only physical values at the be- 
ginning of the Sedov phase. In this paper, SNRs are char- 
acterized by only to, and _Rs- Therefore, our results during 
t > ts do not depend on any assumptions during t < ts in 
our formalism. 

The evolution of the maximum energy of CR electrons, 
Em,c{t), is given by 



-Em,c(t) = min {E 

m , age 



(26) 



Conditions i5m,agc(tc) ~ i5m,cooi(tc) for the free expansion 
phase {t < ts) and £m,cooi(te) = Ein,esc{te) for the Sedov 
phase {ts < t) provide two transition times, tc and te, re- 
spectively (see Figures ([T]) and ((2])). 



SNR at t = te. The maximum energy of runaway CR elec- 
trons is given by 



Eui,c{ic^ 



-Bkno 



En 



3/2 



E,. 



' 5.21 X 10" eV ( for rj^ = r^g.^eo ) 
1.39 X 10^^ eV ( for cx ) (3Q) 
2.66 X 10" eV ( for oc u^i, ) 



where we assume a — 2.6, Equation (|22p for ob and nor- 
malizations of Equation (|25p . Therefore, SNRs can produce 
the Galactic CR electrons up to about 0.3 — 50 TeV which 
depends on the evolution of the magnetic field. Interest- 
ingly, the maximum energy of runaway CR electrons be- 
comes smaller as the maximum energy of CR protons (-Ekncc) 
becomes larger. 



2. 5. 1 The case of tc < ts 



2.5.2 The case of t^ > ts 



We find tc < ts for typical SNRs with a « 2.6. In this 
case, the maximum energy of CR electrons, Bm,e, can be 
expressed by 



Em.c{t) 



where tc is given by 




{t < tc) 
[tcfit^ ts) 

its^t< tc) 
{tc < t) 



(27) 



ts 



1.28 yr 

cc 



^g,frcc 



U015 



5 eV 



ts 



200 yr 



Rs 
2 pc 



(28) 



and te is the start time of escape of CR electrons and given 

by 



te 



ts 



V -C/knee / 



8'7g,frcG')accei5j^jjj,j, 

9.70 X 10^ yr ( for % = %,frcc ) 

3.91 X 10^ yr ( for B^ oc ufi^ ) 

7.39 X 10^ yr ( for B^ oc ) 



(29) 



where we assume a = 2.6, Equation (|22p for as and nor- 
malizations of Equation (|25p . 

In Figure [1] we show the schematic picture of the evo- 
lution of the maximum energy as a function of time (red 
dashed line) for rj^ — 7?g,frco (see also Table 1). In Figure [S] 



show the same figure as Figure[TJ but for B^ 



oc u^h or ii„ 



Initially, the maximum energy increases linearly with time 
(t < tc), then the maximum energy is constant until the be- 
ginning of the Sedov phase (tc < t < ts). During ts < t < tc, 
the maximum energy of CR electrons increases with time 
for rig = ?7g,frcc (Figure [1]) and decreases for B^ oc u^j^ or m^j^ 
(Figure |2]). Finally, CR electrons start to escape from the 



If a is smaller, tc can be larger than te. In this case, Eu 
can be expressed by 



Skncc (t^) 
En..S 

Ekncc{j^) 



(t tc) 

(tc s; t ^ ts) 

(ts t < ts) 

(tB < t < tc) 
(tc < t) 



, (31) 



where we assume a > 1, otherwise i?m,cooi is always smaller 
than -Em, CSC, that is, CR electrons can not escape from the 
SNR until the SNR shock disappears. The start time of es- 
cape of CR electrons, te, for te > tB is obtained by the 
condition i?m,csc(t) = -Em,cooi(t) for ts < t, 



t< 



=tsr 



E„ 



V -Eknc 



(32) 



2.6 Cooling break of CR electrons 



The evolution of the CR spectrum inside S NRs is nec- 
essar y to understand those of runaway CRs (|Ohira et al.l 
I2OIOI ). CR electrons accumulate in SNRs with time, and 
they cool down due to synchrotron radiation during the 
accumulation. Hence, the numbers of cooling CR electrons 
and non-cooling CR electrons are approximately given by 
tcooi,c{E)q{E) oc E~^q{E) and tq{E), respectively, where 
q{E) is an injection spectrum per unit time. Therefore, the 
CR electron spectrum inside SNRs has a broken power law 
form (this is, so called, the cooling break). The break energy, 
i5b(t), is obtained by tcooi,e(£'b, t) = t. From Equations pO|l 
and p8p . the break energy can be represented by 



Eh{t) — -Eb,s 



(t < ts) 

(ts s; t s; tB) , (33) 



where -Eb,s is the break energy at the beginning of the Sedov 
phase (t = ts) and given by 



E, 



b,S 



= £k 



E„ 



-Ekn 
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1.29 X 



(— 



V ts V f Rs 
2 pc 



5.5 eV 



y 200 yr 



(34) 



In Figures [T] and (2] we show the evolution of the break 
energy i5b(i) as a function of time (Cyan line). During 
Eh(t) < Ein,e{t), the spectrum of CR electrons inside the 
SNR has the cooling break. Interestingly, the break energy 
becomes the same as the maximum energy of CR electrons 
(-Em.c = Eh) at the start time of escape of CR electrons 
{t — te). The reason is as follows. Characteristic ener- 
gies, £m,ago(t), Em,cBc{t), Ecooi,c{t), aud Eh{t) are obtained 
from the conditions iacc(-B,t) = t, tcsc{E,t) = facc(-E,t), 
tacc(S,t) = tcooi,e(-B,f), aud tcooi.c{E,t) = t, respectively 
(see also Table [T]). The SNR age is approximately the same 
as the escape time of just escaping CRs {t ~ tasc) because of 
£"111, ago (i) ~ i5m,csc(i) duriug the Sedov phase. CR electrons 
start to escape when the cooling-limited maximum energy 
becomes the same as the escape-limited maximum energy 
(i?cooi,e = Sm.Gsc), that is, the cooing time of CR electrons 
becomes the same as the escape time (tcooi.e ~ icsc). There- 
fore, the SNR age is approximately the same as the cooling 
time {t ~ tcooi,c) at tc. This is the condition to derive the 
break energy. Hence, all the characteristic energies have the 
same energy, i5m,e(ie), at te- 

Namely, the spectrum of CR electrons inside the SNR 
has a cooling break before CR electrons escape {t < tc). 
However, it does not appear while CR electrons escape (t > 

te). 



2.7 Spectrum of runaway CR electrons 

As mentioned above, the spectrum of runaway CR electrons 
is unaffected by synchrotron cooling below En-i{tc)- More- 
over, the spectrum of runaway CR electrons does not depend 
on the magnetic field evolution except for the maximum en- 
ergy of runaway CR electrons. 

According to lOhira et all (|2010h : lOhira fc lokal (|2011^ ■ 
the energy spectrum of runaway CR electrons, dN^sc/dE, 
can be expressed by 



d£; 



AE 



(35) 



where A is the normalization factor, s ~ 2 is the index of 
the energy spectrum of CR electrons inside the SNR and 
/3 is a parameter to describe the evolution of the number 
of CR electrons inside the SNR {dN/dE oc t'^E'' where 
dA''/di5 is the CR spectrum inside the SNR). Lower-energy 
CRs can be produced for a longer time than high-energy 
CRs, so that the time-integrated spectrum becomes softer 
than the instantaneous one. Note that P is not understood 
well for CR electrons and protons because their injection 
processes for particle acceleration have not been understood 
well. In this paper, we assume that CR electrons and protons 
have the same value of /3 = 0.6 during the Sedov phase and 
/3 = 3 during the free expansion phase where we consider the 
therm al leakage model as the injection model jOhira et al.l 
I2010I ). For s = 2.0, a = 2.6 and /3 = 0.6, the spectral index of 
runaway CRs is s + (3/a = 2.23, which is consist ent with the 
index expected as the source of Galactic CRs l|Ohira et al.l 
I2OIOI ). 




10 10 
Energy [ eV ] 

Figure 3. Diffusion length -Rdiff (Equation l|40|l ') and escape ra- 
dius -Rose (Equation II17I I witli a = 2.6) as functions of electron en- 
ergy. The black line shows the escape radius. The red, green, and 
blue lines show the diffusion length at t = 10^, Iff*, and 10^ yr, 
respectively. The solid line shows the case of x = 0.01 and 5 = 0.3 
and the dashed line shows the case of x = 1 and 5 = 0.6 where x 
and 5 are the normalization factor and the energy dependence of 
the diffusion coefficient of the ISM, respectively (Equation (|4TJ). 



3 SPATIAL DISTRIBUTION OF RUNAWAY 
CR ELECTRONS AROUND AN SNR 

In this section, we derive the distribution of runaway 
CR electrons at a given distance r from the SNR cen- 
ter and at a given time t from the explosion time, 
f{t,r,E) [particles cm~'^ eV~^]. CRs propagate around 
SNRs after escaping from the SNRs, so that CR halos are 
produced around the SNRs. Suppose that the system is 
spherically symmetric and the diffusion coefficient in the 
ISM, Dism{E), is spatially uniform. Then, we solve the dif- 
fusion equation given by 



df' 



DisM{E)Af + -^(PiE)f) = q4t,E,r) 



(36) 



where P{E) and qs{t,E,r) are the energy loss rate of CR 
electrons and the source term of CR electrons, respectively. 
Cons i dering the escap e process (|Ptuskin fc Zirakashvilil 
I2OO5I : lOhira et al]|2010l 'l. CRs with an energy of E escape 
from an SNR at t = T'esc(-B). Let dNcsc/dE be the total 
spectrum of runaway CR electrons, that is. 



diVesc 

dE 



dt / d^r qs{t,r,E) = AE 



(37) 



Then, the Green function, which is the solution in the case 
of an instantaneous point source, (ira(t, -E,r) — 5{r)5{t — 
tcsc(E))dNcsc/dE, is lAtovan et al.ll 19951 ) 



fpoint{t,r,E) 



P{Eo) dN,sc{Eo) 



7r3/2i?diff(t,-B)3 P{E) dE 
where Eo{t,E) is the initial energy defined by 



Eo{t,E) 



de 



~ t — TcBc{Eo) , 



(38) 



(39) 



and the diffusion length while CR electrons cool from Eo to 
E is 
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RME) = 2 



Eo(t,E) 



-Pism(';) 

\pi^)\ 



de 



(40) 



10' 



The diffusion coefficient around an SNR has not been 
understood well so far. Affhough there is no guarantee that 
the diffusion coefficient has a power law form, Dism oc , 
(e.g. lFuiita et al ] |2010l .[ 201 ll ) , such a diffusion coefficient can 
explain Kanini a.-ray spectra of rn i ddle-aged SNRs observed 
by Ferrrn (e.g. lOhira et al.ll201ll : \u fc Chenll201lh . There- 
fore, we assume the diffusion coefficient to be 



Djsm{E) = 10^* X 



E 



10 GeV 



(41) 



The values of S and x uncertain. T he Galactic mean 
values are 5 « 0.3 — 0.6 and x ~ 1 (Be rezinskii et ahl 
Il990h . However, x is expected to be much smaller than unity 
(x ^ 0.01) around SNRs because CRs ampl i fy magnetic 
fluctuations a,round SNRs llFuiit a et al.' I2OO9I. I2OI0I. I20T1I : 
iGiuhani et al.ll2010l : [ibrres et al. 2010: Li fc Chen|[201ol ). 

We find that the effect of finite source size is impor- 
tant. To see this, we show, in Figure O the diffusion length 
RdisiE) (color lines) and the escape radius Resc{E) (black 
solid line); the former is obtained from Equations (|16|) . (|39|) 
and (|40|l . while the latter is obtained from Equation (|17|) . 
where we consider synchrotron cooling with B — 3 fiG and 
inverse Compton cooling wit h the Galacti c radi ation field 
provided by the 8 kpc model of lPorter etHI (120081 ). We fully 
conside r the Klein-Nishina effect fo r inverse Compton scat- 
tering (jBlumenthal fc Gouldlll970l '). As shown in Figure E] 
Rdis{E) has low- and high-energy cutoffs. The lower cutoff 
corresponds to the energy of just escaping CR electrons at 
the SNR age, t = tcsc(So). The higher cutoff shows the max- 
imum energy of runaway CR electrons at the SNR age t. It 
is found that for t = 10^ yr and x ~ 0.01 (red solid line), 
CR electrons above 50 TeV escape from the SNR and the 
diffusion length of runaway CR electrons is smaller than the 
escape radius. Similarly, for t = 10* yr and x ~ 0.01 (green 
solid line), CR electrons above ~ 100 GeV escape from the 
SNR and the diffusion length of runaway CR electrons with a 
few TeV is comparable to the escape radius. Therefore, the 
point source approximation is not well, that is, the finite- 
ness of the source size is important. We also remind that 
even if SNRs produce CR electrons with 100 TeV, the SNRs 
should be younger than about lO'* yr (cooling time) and lo- 
cate within about 500 pc from the Earth to observe the CR 
electrons directly. 

Now, taking i nto account the fa ct that CRs escape from 
the SNR surface ijOhira et aLlbOld ). the source term is re- 
placed with 



1 dAfee 

Anr^ dE 



'S{r- R.sc{E))5{t-T,sc{E)) 



(42) 



where Rcsc{E) is the radius where the CRs with an energy 
of E escape fro m the SNR. Then, we find the solution to 
equation ([36]) as l|Qhira et al.ll201ll ) 



/oxt(t,r,£) 



( r-J?.c3c( 

V «ditt(t. 



r + ReBc(-E0) ^ 
«ditf(*.-B) ; 



47r3/27?diff(t,£)i?csc(So)r 

P{Eo) dAesc(£o) 



P{E) 



dE 



(43) 
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Figure 4. Radiation spectra for d = 3 kpc, ?7cR,o = 10 erg, t = 
10'* yr, X = 0.01, and 5 = 0.6 where x and S are the normaliza- 
tion factor and the energy dependence of the diffusion coefficient 
of the ISM, respectively. The red and blue lines show synchrotron 
radiation and inverse Compton scattering, respectively. The thick 
solid Une shows a contribution from runaway CR electrons out- 
side the SNR {Rsh < P ^ Psizc) (Equation l|47| l'). the thin solid 
line shows a contribution from trapped CR electrons which have 
not escaped (Equation JSO}), the long dashed line shows a con- 
tribution from runaway CR electrons which have been caught 
up by the shock (Equation 1148 fi t, and short dashed lines show 
a contribution from runaway CR electrons which have not been 
caught up by the shock but are inside the SNR projected on the 
sky (p ^ Rsh) (Equation 1149 llh respectively. The thick and thin 
black lines show CTA 50 h l ICTA consortiumll2010h and LHAASO 
(where the MAGIC-II- type telesc opes improves the low energy 
LHAASO sensitivity) llCaol I2OI0I) integral sensitivities of point 
sources, respectively. 



The finiteness of the source size is important in the vicin- 
ity of the escape radius (r ~ -Rose) for Rcsc S> -Rdiff- This 
means r — R^sc < -Rdiff a nd r + Resc > -Rd iff, so that /cxt is 
approximately given by (|Ohira et al.ll20l"lh 



/cxt(t,r,S) 



PjEo) dAcsc(go) 



47r3/27?diff(t,-E)-Resc(£;o)r P(E) dE 



.(44) 



4 RADIATION SPECTRA DUE TO RUNAWAY 
CR ELECTRONS FROM SNRS 

In this section, radiation spectra from runaway CR elec- 
trons are calculated. They should be calculated using their 
volume-integrated distribution function. At first, we con- 
sider the energy spectrum of runaway CR electrons con- 
tributing to the radiation spectrum, dNproj/dE. Emission 
far outside the SNR is too dim to be detected because the 
number of emitting electrons is too small. Taking into ac- 
count an observation sensitivity, we should consider the pro- 
jected distance psizc from the SNR center on the sky. We 
adopt Psizc = 1.2-Rsh. Then, we obtain 

dA'^Dioi 



dE 



dp2-Kp dz f,^t{t,x/z^+p'',E) , (45) 



where the z-direction is along the line of sight. In the 
energy region satisfying the condition p <^ Rdiff{t, E), 
the z-integration becomes approximately /point J?diff. Then, 
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Figure 5. Time evolution of radiation spectra from the SNR (0 ^5 
P ^ Psize = l-2iish) for d = 3 kpc,i7cR,o = 10*** erg,x = 0.01 
and 5 = 0.6. The thin and thick color lines show total synchrotron 
spectra and total inverse Compton spectra, respectively. The red, 
green, and blue lines correspond to spectra at the age of 10^, 10*, 
and 1 0^ yr, respectively. The thick and thin black lines show CTA 
50 h ||CTA consortiumll2010l 'l and LHAASO (where the MAGIC- 
Il-type tele scopes improves the low energy LHAASO sensitivity) 
l|Cac I2OI0I ) integral sensitivities of point sources, respectively. 




Energy [ eV ] 

Figure 6. The same as Figure [S] but for 5 = 0.3. 
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Figure 7. The same as Figure [S] but for X = 1- 



dNpioj/dE can be written by 

2 



Psizc 



dE \ RiiH{t,E) J dE 



(46) 



where radiative losses are neglected. 

We further divide runaway CR electrons into three 
groups as in the following: 

• The energy spectrum of runaway CR electrons out- 
side the SNR projected on the sky (_Rsh < p ^ Psizo), 
roj,out/d-E, is obtained by 



dNp 



dp 27rp / dz /e.t {t, ^z-^+p^, E) , (47) 



dE 



• Some of runaway CR electrons are caught up by the 
SNR shock because the SNR shock expands in the halo of 
runaway CR electrons. Note that most runaway CR elec- 
trons are outside the SNR because the diffusion length in- 
creases with time faster than the SNR shock. The energy 
spectrum of runaway CR electrons which have been caught 
up by the shock, dA'caugw/d-B, is obtained by 



diV, 



caught 



dE 



dr 47rr /cxt(t,r, E) 



(48) 



• The energy spectrum of runaway CR electrons which 
have not been caught up by the shock but are inside the SNR 
projected on the sky (p ^ Rsh), dA'proj.in/d-E, is obtained 
by 



diVprc 



dATproj dA^pro 



dTVc. 



dE 



dE 



dE 



dE 



(49) 



Finally, we consider trapped CR electrons. CR electrons 
with an energy smaller than the escape-limited maximum 
energy Em,csc(t) are still trapped in the SNR. Their energy 
spectrum for t > tc, dA^trap/diJ is given by 



diVt. 



dE 



A{Em.csc {t)y 



E 



E„ 



(50) 



In order to calculate each energy spectrum described 
above, several model parameters in section 2 and 3 are 
fixed. Throughout this section, we only consider the model 
of ?7g = Tyg^frcc for the magnetic field evolution, that is, the 
start time of escape of CR electrons is te = 9.70 x 10 yr. 
In our formulation shown in Section [2l the spectral index 
of runaway CR electrons does not depend on models of 
the magnetic field evolution, while the maximum energy 
of runaway CR electrons, -Em,e, and the start time of es- 
cape of CR electrons, to, do depend. We adopt Rs = 2 pc, 
ts = 200 yr, Bism = 3 pG, rj^cc = 10, s = 2.0, a = 2.6, 
/3 = 0.6, d — S kpc, nisM = 1 cm~^, and UcR,e = 10** erg, 
where d, nisM, and UcR,e are the source distance, the num- 
ber density of the ISM, and the total energy of runaway 
CR electrons from 1 GeV to -Bknee, respectively. We calcu- 
late inverse Compton scattering wit h the Galactic radi ation 
field provided by the 8 kpc model of lPorter et all (|2008l ) and 
synchrotron radiation. We fully consi der the Klein-Nishina 
effect for inverse Compton scattering (iBlumenthal fc Gouldl 
Il970h . We use the downstream magnetic field Bd = 4B for 
synchrotron radiation emitted inside the SNR. 
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For given distributions of CR electrons (Equations (|47[l - 
H50[)'). radiation spectra are calculated. Figure [4] shows the 
spectral components of four CR electron groups given by 
Equations (gTll-lO, where t = 10* yr (> U), x = 0.01, and 
5 = 0.6 (x and S are the normalization factor and the energy 
dependence of the diffusion coefficient: see Equation (|4ip '). 
We also show in the figure the integral sensitivities of point 
sources of the future Te V gamma-ray tel escopes, CT A 50 h 
(|CTA consortium! |2010| ) and LHAASO (|Cad |2010| ) (black 
lines). The thin and thick lines show the radiation spectra 
coming from the interior of the SNR region projected onto 
the sky (that is p ^ Rsh) and from the exterior of the SNR 
(i?sh < p ^ psizc — 1.2i?sh) (Equation (|47)) '). respectively. 
The thin-solid lines show spectra emitted by trapped CRs 
which have not escaped from the SNR yet (Equation (|50p ). 
The long and short-dashed lines show the spectra emitted 
by runaway CR electrons which have been caught up by the 
shock (Equation (|48|) ') and have not been caught up by the 
shock (Equation (|49|l '). respectively. 

One can see from Figure [4] that maximum-energy pho- 
tons of synchrotron radiation are emitted by runaway CR 
electrons which have been caught up by the shock. This is 
completely different from the standard picture that maxi- 
mum energy synchrotron photons are produced by CR elec- 
trons accelerating at the SNR shock. Note that this new 
picture can be applied after CR e lectrons start to escape 
(t > te). In previous works (e .g. Berezhko fc Volkl |2010| : 
IZirakashvili fc AharonianI l2010l : lEdmon et al] l201ll '). syn- 
chrotron emission comes from CR electrons advected down- 
stream at any time, that is, there is no halo of runaway CRs. 
In our model, the diffusion coefficient rapidly increases with 
time compared with the previous work. Therefore, in our 
model, CR electrons can escape from inside SNRs over ad- 
vection and the halo of runaway CRs is produced. 

Radiation spectra from runaway CRs are softer than 
that from trapped CRs. This is because higher-energy CRs 
are more diluted by diffusion and the spectrum of runaway 
CRs becomes softer than that of trapped CRs (see Equa- 
tion H35|) ). From Equations (|38p and (|48|l . the energy spec- 
trum of caught-up CR electrons can be written by 



caught 



dE 



f -Rsh \ 



cx 
oc 



-3.13 



'dN^ 
dE 

s + ^ + 



(51) 



where we assume s = 2.0, /3 = 0.6, a = 2.6 and 5 = 2.6 at the 
last equation. Therefore, the energy spectrum of synchrotron 
radiation from caught-up CR electrons becomes flat (long- 
dashed line). In addition, we predict that inverse Compton 
scattering outside the middle-aged SNR (thick blue line) will 
be observed by CTA and LHAASO if the SNR produces 
sufBcient CR electrons (?7cr,g = 10** erg) and x = 0.01. 

In Figures [SH71 we show time evolutions of radiation 
spectra around the SNR (p < psizc = 1.2Rsh) for three pa- 
rameter sets of the diffusion coefficient of the ISM, (x, 5) ~ 
(0.01,0.6), (0.01,0.3), and (1,0.6). Integral sensitivities of 
point sources of th e future ground-based Cherenkov tele- 
scopes, CT A 50 h ||CTA consortium! I2OI0I ) and LHAASO 
l|Cad l20ld l (black lines), are also shown. The thick solid 
and thin solid lines show total synchrotron radiation and to- 
tal inverse Compton scattering, respectively. Figures [5] and 




10" 10 10 
Energy [ eV ] 

Figure 8. Time evolution of gamma-ray spectra from the SNR 
(with the projected distance Psizc < 1.2Rsii) for C/cR,p = 
10^" erg, nisM = 1 cm~^,x = 0.01,5 = 0.3, and other param- 
eters are the same as Figure [5] The thick and thin lines show the 
total inverse Compton spectra and the gamma rays {n^ — >■ 27) 
from runaway CR protons, respectively. The red, green, and blue 
lines correspond to spectra at the age of 10^,10*, and 10^ yr, 
respectively. The thick and thin black lines sho w CTA 50 h 
||CTA consortium! I2OI0I ) and LHAASO l!Cad!2010h sensitivities, 
respectively. 



[6l (small diffusion coefficient case) show that even though 
the SNR age is 10^ yr, the future Cherenkov telescopes can 
detect radiation from runaway CR electrons. 

Figures [SI and [71 show that radiation from runaway CR 
electrons depends on the normalization of the diffusion co- 
efficient, X- As shown in Equation (|46|l . the flux of in- 
verse Compton scattering is large for small x- Therefore, 
if X < 0.1, radiation due to runaway CRs from middle- 
aged and old SNRs can be potentially observed by CTA and 
LHAASO. Similarly, as shown in Equation (|5H) . the flux of 
synchrotron radiation from caught-up CR electrons is large 
for small x- It should be noted that if the evolution of the 
magnetic field is oc ufj, or u^jj, the maximum energy of 
caught-up CR electrons is too small to produce synchrotron 
X-ray photons (Equation (|30|l '). However, if the magnetic 
field at the SNR sho ck is suddenly am plified by interactions 
of molecular clouds (|lnoue et al.!!2009l ) , X-ray photons could 
be produced by synchrotron radiation from caught-up CR 
electrons. 



5 DISCUSSION 

In this paper, taking into account escape of CR electrons 
from SNRs, we have investigated the evolution of the maxi- 
mum energy of CR electrons, as well as that of the magnetic 
field and the maximum energy of CR protons. We found the 
foUowings: 

• CR electrons start to escape from SNRs later than CR 
nuclei (Equation P^ l. 

• SNRs can produce CR electrons up to about 0.3 — 
50 TeV, which depends on the evolution of the magnetic 
field (Equation 

• If the magnetic fleld decays slowly, the maximum energy 
of runaway CR electrons becomes small (Equation (|30|) ). 
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• The maximum energy of CR electrons is anticorrelated 
with that of CR nuclei (Equation (|30|l '). 

• The spectrum of runaway CR electrons is unaffected by 
the cooling break of synchrotron radiation (Section H2.6[) V 

In addition, we calculated the spatial distribution of run- 
away CR electrons around SNRs. Using this distribution, 
we calculated synchrotron radiation and inverse Compton 
scattering from the runaway CR electrons around SNRs. We 
found the foUowings: 

• Inverse Compton radiation of the runaway CR elec- 
trons can be observed by future ground-based Cherenkov 
telescopes, such as CTA and LHAASO. Because the cos- 
mic microwave background is inevitably present, CTA and 
LHAASO can impose tight constraint on the amount of CR 
electrons. 

• Maximum-energy photons of synchrotron radiation are 
emitted by runaway CR electrons which have been caught 
up by the shock. 

• Radiation spectra from runaway CR electrons depend 
on the diffusion coefficient of ISM. 

Although there are many parameters in determining the 
evolution of the maximum energy (_Rs, ts, -Bism, -Eknoo, '7acc, 
'7esc, ??g,tree, Q and qb ; see Table 1), all except qb have typi- 
cal or appropriate values. The value of qb , which is a param- 
eter to describe the evolution of the magnetic field, is highly 
uncertain. We have derived some relations for the general 
Qb and especially discussed three cases. We have assumed 
that the end of the Sedov phase is 200 yr x 10 ^'^ ~ 63 kyr. 
On th e other hand, iTruelove fc McKeel (| 19991 ') and iPetrukI 
l|2005ll thought that it is 12 kyr and 30 kyr, respectively. 
Even when we take these values, 12 or 30 kyr, our results 
do not change significantly. In this paper, we have consid- 
ered only one fiducial model of SNRs (ts = 200 yr and 
Rs = 2 pc). It is also an interesting future work to ex- 
amine the individual SNR detected by H.E.S.S., MAGIC, 
and VERITAS. It should be noted that the age of most 
young SNRs is about the beginning of the Sedov phase and 
smaller than the start time of es cape of CR electrons, tha t 
is, t r^ts <U (see the table 1 of lTruelove fc McKeelll999t ). 
Therefore, CR electrons have not escaped from those SNRs. 

We consider three evolution models of the magnetic 
field in Section [2l Different models predict different evolu- 
tions of the maximum energy limited by synchrotron cool- 
ing during to, < t < min{tc,tB}. Maximum photon ener- 
gies of synchrotron radiation and inverse Compton scatter- 
ing depend on evolution models of the magnetic field. Fur- 
thermore, the maximum energy of runaway CR electrons 
also depends on the evolution model of the magnetic field 
(0.3 - 50 TeV). Figures [5] H and [7] show that CTA and 
LHAASO can observe the maximum value. Hence, CTA and 
LHAASO can provide constraints on evolution models of 
the magnetic field. To make a more accurate radiation spec- 
trum, we should consider escaping and trapped CR spec- 
tra in detail (Caprioli et a l. 2009; Rcvillc ct al. 20o3) and 
solve SNR dynamics and t he diffusion-convection equation 
llZirakashvili fc Aharonian|[201Q '). 

Fermi and H.E.S.S. show the spectral break or the cut 
off of the CR electron spectrum at around 1 TeV. The mag- 
netic field evolution model of oc could explain the 
break as the maximum energy of CR electrons that can 



escape from SNRs (see Equati on (1301) ). This model is sug- 
gested by other observations (iVinla2008l ) and by JSell (200J). 
A local source may contribute to high energy CR electrons 
and the maximum energy of CR electron s may be limited 
by c o oling during the propagatio n (e.g. iKobavashi et al.l 
12004 iThoudam fc Horandell [20121 ). Even so, we can rule 
out any models that predict a smaller maximum energy 
of CR electrons than that of observed CR electrons. A 
more tight constraint on the magnetic field evolution will 
be given by future observations of the CR electron spec- 
trum, AMS-02 and CALET. To make a more realistic pre- 
dict ionj_we_^h£uIdpCT^ time-dependent fluid simula- 
t ion s JSchure et aT bold) and consider the CR back reaction 
fe.g.'Prurv et al.'ll989l:lEllison et al.l ll990l : lBerezhko fc Volkl 
[1997; Blasi 2002; Kang ct al. 200^ 

In this paper, we have focused on escape of accelerated 
electrons. SNRs produce CR protons as well as CR electrons, 
so that the hadronic gamma-ray radiation due to the neutral 
pion decay should be compared to the leptonic gamma rays. 
In figure[8l we show an example of the gamma-ray spectrum. 
We adopt Ucn,p = 10^° erg,nisM = 1 cm~^,x = 0.01,5 = 
0.3, and otherparameters are the sa me as in the previou s 
section. We use the code provided bv iKamae et al.l l|200d ): 
iKarlsson fc Kamad (|2008 ') Q 

to obtain gamma-ray spectra 
from CR protons. Note that if the ratio of CR electrons to 
CR protons is (7cr,c/(7cr,p = 10"^ and the number density 
of the ISM is nisM ~ 1 cm^'^, inverse Compton scattering 
due to runaway CR electrons dominates over the hadronic 
component at the TeV energy region. The maximum pho- 
ton energy and its flux decrease with time for synchrotron 
radiation and inverse Compton scattering because of cool- 
ing and diffusion. On the other hand, the maximum energy 
of gamma rays from CR protons does not change because 
the cooling time of CR protons is longer than the SNR age. 
Therefore, radiation from CR protons can dominate over 
that from CR electrons at the 100 TeV energy region for 
middle-aged and/or old SNRs. 

We comment on a possible origin of unidentified very- 
high-energy gamma-r ay sources (so called, TeV-unlDs) dis - 
covered by H.E.S.S. l| Aharonian et all |2005| . l2006bl . |2008| ). 
Figures [Hlini and [7] tell us that old SNRs {t ~ lO'^ yr) could 
be the origin of TeV-unlDs. These figures show that inverse 
Compton scattering from runaway CR electrons can be ob- 
served in the TeV region, but synchrotron radiation can not 
produce X-ray photons for old SNRs. The maximum energy 
of runaway CR electrons decreases with time owing to syn- 
chrotron and inverse Compton cooling and the m aximum 
energy is a few TeV at t — 1 0^ yr (see also lYamazaki et al.l 
12009 : lloka fc Meszarod [2OI0I ) . Electrons with energies of a 
few TeV can produce TeV gamma-rays by inverse Compton 
scattering with optical photons, but can not produce X-ray 
photons by synchrotron radiation with the magnetic field 
of 3 /^G. This is the same reason why pulsar win d nebulae 
(PWNe) are plausible candidates for TeV-unlDs l|de Jagerl 
l2008l ). Radiation from runaway CR electrons depends on the 



^ Note that in Figure [S] the softening around 30 TeV of gamma 
ra ys from CR protons i s due to the limitation of th e code provided 
bv lKamae et al.l ll20od) : lKarIsson fc Kamaj 1I2OO8I ). The code cal- 
culates gamma-ray spectra from CR protons with energies up to 
512 TeV. 
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number of CR electrons, the difFusion coefficient around the 
SNR, the SNR age, and the source distance. For example, if 
the value of x is larger than 0.1, the flux of inverse Compton 
scattering from runaway CR electrons around an SNR older 
than 10* yrs with the distance 3 kpc is smaller than the 
sensitivity of H.E.S.S., so that the number of currently de- 
tected TeV unlDs is not so large. However, even in this case, 
we expect that future observations by CTA will increase the 
number. 

In addit ion, we comment on e scape of CR electrons 
from PWNe. iKawanaka et al.l (|201ll ) investigated CR elec- 
tron spectra from a young pulsar embedded in the SNR, but 
did not take into account synchrotron cooling while CR elec- 
trons cross the SNR shell. If the diffusion escape ti me from 
the SN R is longer than the cooling time, Kawanaka et al.l 
l|201ll ) is not valid. However, CR electrons are not pro- 
duced at SNR shocks for PWNe. That is, the maximum 
energy of CR electrons produced in PWNe could be larger 
than that produced at SNR shocks. In this case, CR elec- 
trons produced in PWNe can escape from SNRs as long 
as the diffusion escape time is smaller than the cooling 
time. This depends on the evolution of the magnetic field. 
IKawanaka et al.l l|201ll ) proposed that the low energy cutoff 
of CR electrons escaping from nearby sour ces could be ob- 
served in the above TeV band. Therefore, iKawanaka et al.l 
l|201ll ) is complementary to this paper that investigated the 
high energy cutoff. 

Finally, we emphasize that escape of accelerated par- 
ticles is common in other high energy particle accelerators 
(e.g. microquasars, gamma ray bursts and active galactic 
nuclei), so that radiation of runaway accelerated particles 
from the accelerators can also be expected. 
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